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ABSTRACT 

The MPMS mechanism possess two revolute degrees-of-freedom and allows the 
user to measure the mass, center of gravity, and the inertia tensor of an unknow 
mass. This paper develops the dynamics of the Mass Properties Measuremen 
System (MPMS) from the Lagrangian approach to illustrate the dependency 
the motion on the unknown parameters. 


93 



Kennedy Space Center 
MPMS Dynamics 


Aug 6, 1993 


Keith L. Doty 
NASA Faculty Fellow 


4 


1. INTRODUCTION 

The Mass Properties Measurement System (MPMS), illustrated in Figure 1.1, 
consists of a serial kinematic mechanism with two intersecting revolute axes, z l 
and z 2 , that intersect with fixed angle a v The joint angles and rates denoted by 

Gj, 0 2 and 0 l5 0 2 , respectively, turn about the respective joint axes. 



Figure LI Mass Properties Measurement System (MPMS). 


The function of the MPMS is to measure the mass, center-of-gravity and second- 
order mass moments of an unknown mass placed on the turntable. The mass and 
center-of-gravity can be determined from static measurements. The inertia tensor 
must be determined from the dynamics. 

The structure of the Mass Properties Measurement System (MPMS) lends itself to 
a straightforward dynamics analysis using the Lagrangian approach. The 
dynamics analysis assumes ocj as a parameter. 

The symbolic software package Mathematica was used to produced and verify 
the equations presented here. In the analysis to follow, the standard Denevit- 
Hartenberg kinematic parameters of the MPMS were employed. 
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2. MPMS KINEMATICS 

Table 2.1 lists the Denevit-Hartenberg kinematic parameters for the MPMS 
mechanism. 

Table 2.1 : Kinematic Parameters for the MPMS Mechanism 


Joint 

d 

0 

a 

a 

1 r 

0 

0i 

0 

OCi 

2 7l 

0 

02 

0 

0° 


From the DH-parameters of the MPMS mechanism listed in Table 2.1, the two 
link transforms compute to 

/ 0, 


Li = 


Ri 


0 3 


0 0 0 1 


L 2 = 


f x Ro o ^ 


0 0 0 1 


(2- la) 


f 

ci 

- si 

a i si \ 


f c 2 

-S2 

0 

\ 



Sl 

Xi C! 

-0\ ci 

l R 2 = 

S2 

c 2 

0 


(2- lb) 


0 

°1 

*1 ) 


V 0 

0 

1 

) 



where °Rj = 


and c, := cos(0j) , Sj := sin^), X; := cos(Oi) and O; := sin(«i). The forward 
kinematics transform of the MPMS equals 

°T, = Lj L 2 (2 ' 2) 


which computes to 


0, 


°T ? = 


0 


r 2 

0 0 0 1 


(2-3) 


where 
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j r 2 = 


c 1 c 2" t 1 s 1 s 2 - c 1 s 2‘ t 1 s 1 c 2 a l S 1 


a l s 2 


° 1 c 2 


(2-4) 


The forward analysis presented here serves as reference. The rotation part 
indicates how to change frame F 2 vector representations into frame 1 0 
representations. The dynamics analysis presented later will make use of the 
forward kinematics °T 2 • 

3. MPMS End Frame Jacobian 


The Jacobian of the MPMS relates the joint-rates q 
velocity V = [v T co T ] T of the end-frame, 

V = Jq. 

The Jacobian of the MPMS computes [1] [2] to 

0 0\ 
0 0 
0 0 

2 , ° J 2,2 = <*1 s 2 0 

Ci c 2 0 

v x i 1 y 


0 2 ] t to the frame- 


(3-1) 


(3-2) 


The leading superscript 2 means that this Jacobian is expressed in frame 2 
while the 0 indicates the motion the end-frame, designated by the second 
subscript, is relative to the base frame F 0 of the MPMS. The first subscript 
indicates the frame origin at which the linear velocity is measured. 


For convenience, we write 


0 0 and 2,0 J. 


CiS 2 0 
CiC 2 o 
x, 1 
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4. MPMS DYNAMICS 

The Lagrangian approach to dynamics of the mechanism requires the calculation 
of the kinetic and potential energy of the moving masses. 

Potential Energy Terms 

Assume the center-of-mass vectors for the first and second links equal r,* and 
r 2 *. These vectors are constants in their own frames, 


l ri = 


AxA 

r y i 

\ r zl ) 


and 


2| *2 = 


r y2 

V r z 2) 


(4-1) 


The gravitational field vector in Figure 1.1 equals 

ro \ 


0 g= gc 


°1 

V X 1 ) 


(4-2) 


and the potential energy of a mass m at position p equals 

P = -mg x p 

hence, the potential energies P x and P 2 of the first and second link equal 


T * 0 % Od 1 w* 

Pj = -nq g T r 1 = -mj g Ri *T 


and 


x * O—T 0i> 2_ * 

P 2 = -m 2 g r 2 = -m 2 g K 2 r 2 . 


(4-3) 


(4-4a) 


(4-4b) 


The torques 


associated with changes in the potential energy of the mechanism 


links equal x pel - ^ 
compute to 


3P, 3P 2 

+ a0| and 


dP, dP 2 3P 2 

„ = -L- + -,J = . These torques 

pe2 00 2 d0 2 
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x pel = m l C 1 + V S 1 + a l S 1 } + m 2 a l Sc( r x2 ( C 1 c 2 - X 1 S 1 s 2> 

+ r y2 (" x l c 2 S 1 " C 1 S 2^ + r z2 g 1 S 1 ) (4-4a) 

T pe2 = m 2 G 1 g C { r x2(- c 2 X 1 + X 1 C 1 c 2 - S 1 s 2> + r y2 ( X 1 s 2 “ c 2 S 1 " X 1 C 1 s 2> ) ( 4 " 4b > 

Equations (4-4) will provide the means for measuring the mass and center-of- 
gravity for each link. 

Measuring MPMS Link Mass and Link Center-of-Gravity 

By measuring the balancing torque on joint one and setting the joint angles at 
different angles one can determine the required information about the first-order 
mass moments of each link. This information will allow us to subtract the gravity 
torque terms due to the mass of the MPMS when we desire to measure the mass 
properties of the unknown mass. 

Experiment 4.1 = 0, 0 2 = 0, 


X xl := m l g c r xl + m 2 <3 X g c r x2 

Experiment 4.2 0j = 0, 0 2 = rc/2, 

X x2 G 1 8c r xl — m 2 G 1 Sc r y2 

From these two measurements one computes 

X xl “ X x2 . X xl + X x2 

m 2 f y2 =T^17 and m 2^i = 2<Tlgc 
Experiment 4.3 0j = ft/2, 0 2 = 0, 

X x3 := m l G 1 8c( r yl X 1 + r zl G l) + m 2 G 1 g c (- r y2 X 1 + r z2 G l) 

Experiment 4.4 0j = 7t/2, 0 2 = n/2, 

X x4 := m l G 1 g c ( r yl X 1 + r zl G l) + m 2 G 1 gc(- X 1 r x2 + r z2 G l> 


(4-5) 

(4-6) 

(4-7) 

(4-8) 

(4-9) 
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From Experiments 3 and 4 one can find 


X x3 - T x4 _ 

m 2 r x2 - Gl g c 2 y 2 


(4-10) 


Relation (4-10) does not apply when x, = 0, i.e., when a! n/2. 

From (4-7) and (4-10) we can develop a linear relation in the other unknowns, 

x pel := (-^ - m l r *l <=1 - m 2 r %2 < C 1 c 2 - S 1 

+ m 2 r >2 (-xic 2 s 1 -c 1 s 2 ) ) = m, (r y , x,+ r zl a, ) + m 2 r^ o, (4-11) 

Since the unknowns in (4-11) have fixed coefficients, those unknowns cannot be 
resolved further. 

r (A i\ (A im and (4-11) provide the necessary information for 

terms, k actual mass values and center-of- 

gravity terms do not need to be determined completely. 

MPMS Kinetic Energy Terms 

The total kinetic energy K of the MPMS motion equals K = K, + ^ 2 .^^ 
i = 1 2, equals the kinetic energy of link L, defined with respect to the origin 
define in Figure 1.1. The joint torques x kel and x ke2 required to generate the 
motion, assuming a conservative system, equal, accordmg to Lagrange, 


T kei = a 


9K' 


89; 


8K 

80 i ’ 


i = 1,2. 


(4-12) 


The torque terms associated with each kinetic energy 

"8K; 


d 

T keji _ dtl 


U e i 


8Kj 
89i ’ 


j = 1,2; i = 1,2. 


(4-13) 


sum to produce the total torque for that joint, 

'Ckei = x keli + x ke2i 


(4-14) 
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Since the origin does not translate, the kinetic energy for both links of the MPMS 
equals the rotational kinetic energy, 


Kj = \ ©i <*> Ii ©i, i = 1>2, 


(4-15) 


where co* equals the angular velocity and Ij equals the standard inertia matrix for 
link Lj . In frame F { the matrix I* has constant terms. 


From (3-1) and (3-3), 


i’V = 0! 


(0\ 

v T v 


2 -°co 2 = 2 '°J 


to,2 


A/ 


The kinetic energy of the first link 

K 1= ! 1 »°o ) 1 <-> 1 I 1 1 ’°<o i 

does not depend on either of the joint variables nor 0 2 , hence. 


T kell : dj 


0K, 


00 , 


0Kl .. ,T Ij 

-007 =0 i (O °i Tl) 11 


0 \ 

<?1 

V Tl > 


x kel 2 : “ d 3 


0K, 


00 - 


0Kj 

00o 


= 0 


The kinetic energy of the second link equals 


Ko = 


2~2 


V 

V 0 2, 


2,0 ¥ X 2 t 2,0 1 

<*> J W 2 *2 Jco,2 


V 

V 0 2y 


(4-16) 


(4-17) 


(4-18) 


(4-19) 


(4-20) 


and does not depend upon 0 4 . Hence, the torque contribution by joint one to the 
motion of the second link equals 


Kennedy Space Center Aug 6, 1993 

MPMS Dynamics 


Keith L. Doty 1 1 

NASA Faculty Fellow 


d 

raK 2 i 

3K 2 

d 

raK 2 -i 

x ke21 :— dt 

* 

1 

30 j 

" dt 

_30i. 


2,0 / 2 t 2,0 1 rh 

J(02 *2 


<•> 


V 

• • 

0. 

V *■) 


2,0 / 2| 2,0 1 <•> 

J (02 * 2 J *>’ 2 VOj < 




1 

\Aj 


2,0/ 2t 2,0 
+ J C02 A 2 


'La (i) 


'e ' 


<•> 


e 




(4-21) 


Similarly, the torque contribution by joint two to the motion of link two equals 


x ke22 : ~ dj 


dK 2 


30 


2 J 


3K; 

ae 2 


qxv 2 2,0/ 2 t 2,0 1 

” " J (02 2 ' 


/0\ 

/*• \ 
9l l 

1 <•> 

•• 

w 

09 

2 ) 


(i) <-> e,] + 2 ° J “ 2 ^ 2 °'* w ' 2 

V 2 / 


A)\ 

A) 

1 1 1 <*> 

. 

w 





( 0 2 


a 2 ’°j 


co2 


<#>1 39 2 

From (3-3) compute 


2 ’X,2 = °1 ^ 


^2.0 ¥ 

2t 2,0 1 + 2,0j X 2j j 

*2 J o>,2 + J to2 2 302 


V 

A, 


(4-22) 


0) 
-s 2 0 
0 0; 


and 


a 2>0 j 


ae, - 0l 


(c 2 0\ 

-s 2 0 

U °J 


(4-23) 


and deduce 


2 ^; 2 % 2 ' o j w ,(?)= 2 ' 0j :2 2i 2 2 'X/i)= ( 0)’ 


0\ /^l( s 2 X xz2 + c 2 I yz2^ 


-4% 2 -x 2 (i)=! r 


^zz2 


(4-24) 

(4-25) 
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Relations (4-24) and (4-25) allows us to simplify (4-22), and, coupled with (4-19) 
yields the joint- two, kinetic energy derived torque x ke2 . 


x ke2 


. _ 2 , 0 2 t 2,0 ¥ (^3 

= x ke22 + x kel2 “ J (o2 * 2 J ®’ 2 UJ 


<•> 


V 

e 2 




v e 2. 


9 2 ’°J ^2,0 t 

_ . 2 2,0 1 . 2,0 y X 2x a 

<•> < ae 9 12 V 2 + J ©2 12 ae 2 1 


V 

(®2 


(4-26) 


The joint-one, kinetic energy derived torque T kel equals the sum of (4-18) and 
(4-21), 


2,0 y 2y 2,0, 


x kel = x kel 1 + x ke21 “ ^Tl + ’ ^©2 *2 


fl\ 

% 

n <*> 


{OJ 

d 2 

l ) 


, , 2,0 2t 2,0, , 2,0, v 2, 

+ { *2 J(o,2 + J co2 12 


2,0 y T 2y 2,C 


/T\ 


A <*> 

• * 

W 

02 

2/ 


(4-27) 


where 


Ixl (0 


x l) X % 


(0 \ 

<*1 

V x iJ 


cannot be resolved further by any joint measurements 
Using Mathematica to expand (4-27) and (4-26) yields 


(4-28) 
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T tel = §i I Tl + I„ Tl (9 2 + 01 t t ) + Oil I yy C 2 ffj (0, C 2 - 2 0, 0 2 . 2 ) 
+ 1^ Gi s 2 (2 0! 0 2 c 2 + 01 S 2 ) 

+ lyz ("®2 c 2 - 2 0i c 2 X! + 0 2 2 s 2 + 2 0J 0 2 Xj s 2 ) 

+ I xz (-0 2 2 c 2 - 2 0! 0 2 c 2 Tj - 0 2 s 2 - 2 0J s 2 ) 

- I^y Gi (2 0 X 0 2 C 202 + 0! s 202 ) } 


and 


^ke2 ^zz (®2 


* 

+ 0 1 X 1 )+ Oil I Iy 01 2 c 292 Ol+Ixz(01 2c 2^- 0 1 s 2> 


+ I y2 (- 01 C 2 - 0| 2 X| S 2 ) - l I x , 01 2 O, S 202 + 2 I yy 01 2 CTl S 202 ) ( 4 ‘ 3 °) 

Equations (4-29) and (4-30) yield the essential relations for determining the 
unknown inertia matrix I 2 and the term I T1 from torque, velocity and 

acceleration measurements. 


5. DYNAMICS OF THE UNKNOWN MASS 

The functional objective of the MPMS is to measure the unknown mass properties 
of an object placed on the table of link two (Figure 1.1). The unknown mass is 
rigidly attached to link two during the measurements. Conceptually, this makes 
the unknown object a part of the second link, hence, the combined inertia tensoi 

I equals 


\ = h + *u 


(5-1) 


where 2 I U is the inertia matrix of the unknown mass described about the common 
origin of frame F 0 and Fj . The dynamics analysis in Section 4 applies to this 

problem with I 2 replaced by Thus, the torques derived from the kinetic 
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energy of the unknown mass alone have the form of T kel - ©i I T i in (4-29) and 
x ke2 in (4-30) with the inertia terms replaced by those of the unknown mass, 


T keul : = Izzu *1 (0 2 + T l) + <*1 { !yyu c 2 <*1 (®1 c 2 ‘ 2 0 2 S 2 ) 


• • • • 

+ *xxu a l s 2 (2 01 02 c 2 + 01 s 2) 

+ ^yzu (" 02 C 2 " ^ 01 c 2 T 1 + ®2 2 s 2 + ^ 01 02 T 1 s 2) 


+ ^xzu ("02 2 c 2 ' ^ ©2 c 2 T 1 " 02 s 2 " 2 01 X 1 S 2 ) 


• • « « 

“ ^yu a l (2 ©1 0 2 C 2 0 2 + 01 S 2 02) } 


(5-2) 


and 

*ke2 =I zzu(02 + 01 X l) + <7i{I xyu 0l 2 C 202 ai+ I XZ u(01 2 C 2 'Cl- ©1 S 2 ) 

+ Iyzu (-01 c 2 - 01 2 X 1 S 2> - 2 ^xu 01 2 G 1 ^02 + 2 ^yu ®1 2 a l S 202 ) ( 5_3 ) 

Equations (5-2) and (5-3) yield the essential relations for determining the 
unknown inertia matrix I u from torque, velocity and acceleration measurements. 

6. Measuring the Unknown Mass 

Since the torque equations are essentially the same for the inertia terms of link 
two of the MPMS and the unknown mass, only the experimental technique for the 
latter will be developed here. 

The first question to resolve: 

Will measurements of the torque, angular position, velocity and acceleration of 
joint one provide sufficient data to compute all the inertia parameters of the 
unknown mass? 
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The answer will to this question is No. To verify this claim note ^t potenhally, 
all six independent parameters in I u might be determined from ( 5 - 2 ) given that 

Experiment 6.1: 0 2 = 0 , 0 2 = 0, 

T keol := 01 (I«u V + a l 2 l c 2 2 + ^ s 2 2 ' ^ 5262 1 

- 2 0 1 t [ { I yzu c 2 lizu s 2 ) ^ 
which, in inner product form, equals, 


-W = 9i k <•> 'v (6_1) 

where k is the coefficient vector and I v a six-vector of the independent 
parameters in I u : 


k := 


Tl 2 A 



2 2 
<?1 c 2 


fyyu 

,2 2 


I 

<*1 s 2 
2 

and I v := 

XXU 

I „„ 

'<*1 s 202 


xyu 
^ yzu 

vw 

- 2 X! C 2 

V 2 01 *1*2 J 



( 6 - 2 ) 


If a vector can be found orthogonal to the coefficient vector for all choices of 
MPMS configurations, then, six independent coefficient vectors k cannot be 
“d and therefore, six independent equations in the six unknowns ,n I v 

cannot be obtained. The vector x := [a, 2 -x, 2 V 0 0 0f*0i.«»ch. 

vector- k <•> x = 0 for all possible k . This implies no non-singular 
measurement matrix M can be formed from the set of possible k that will yield 
six independent equations. Since M x = 0 for any M whose rows are constructe 
from different k's, dte matrix M possesses a non-zero vector in ,te nuU : space ^and 
therefore, must be singular. Note that this observation holds for any choice o 
angle 04 between the two joint axes of the mechanism. 

Observe that x cannot be solution to any inertia problem because it forces at least 
one of the diagonal terms of l u negative, an impossible situation physically. 
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Substitute x = I v , this is not a physically possible I v , into the general expression 
for T teol and find that T ke „, = I zzu C! t, 0 2 which proves that, minimally, one 
must measure 02 on the second joint. 

♦ • • 

Experiment 6.2 0j = 0 , 0j =0, 

'W : = I zzu T 1 ^2 + a l I yzu (- ^2 C 2 + *2 S 2 > + a l *xzu ^2 C 2 ' ®2 S 2> (6 ' 3) 

hence. 


^keul = k <#> *8 


(6-4) 


where k is the coefficient vector and I s a three-vector . 

f - \ 


k := 


Tl 0; 


<3 i (“02 ^2 ®2^ ^2 ) 

\^CTi (-02 ^ c 2 ” ^2 S 2) J 


and Io := 


rtzzu^ 

^yzu 

VW 


(6-5) 


By measuring the torque on joint one and the angular position, velocity and 
acceleration on joint two, one can measure three independent equations (6-4) to 

_ n t 2 • n t •• n_t 

compute 1^ I yzu , I xzu . For example, if 0 2 = 4 , 0 2 - 2 » w 2 “ 2 > 

0 < t < a J~2 , and a ^ = tan _1 [V 2 j, the measurement values of (6-3) for t=0, 
t = 1, t = V 2 lead to the measurement matrix 



^ k T (t=0) ^ 


f 0.9069 0 -1.28255 'j 

M = 

k T (t=l) 

(k T (t=%/2 )J 


0.9069 -2.33145 0.517655 
0.9069 -1.28255 4.02925 J 
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whose determinant equals det[M] = -9.13734, proving that M non-singular and 

ft A . . . 

be solved for I s . With I s known, the remaining inertia 


that MI. = 




can 


parameters can 


be determined from the measurements in Experiment 6.1. 


on 


Conclusion ■ If one measures torques and angular « cc f. lerati ° n .J" 

Sr s r iswsaars vs 

inertia parameters of the unknown mass. 


7. Proposed Design for MPMS 

Based upon the analysis in this paper, we note that 

1 . All the mass properties of the MPMS and the unknown mass itself ca^ 
be determined by measurements performed only on joint axis one 

quantities. 

Discussions with Kedron Wolcott on the construction of the current prototype 
indicates 

2. Construction of the bearings for the MPMS for a twist 0 < <*i < 90° 
present significant cost and design penalties. 

A talk with Richard Bennett indicates that 

3. The current design occupies a large volume due to the twist «i = 54° 
between the two axes. 

4. The current design is not easily scalable to handle larger loads. 

Based upon the four considerations above, the a»*or Proposes ’ *^1 the two joint 
axes of the MPMS be oriented at right angles, a, = 90 , as shown in g 
This (lesion eliminates the problems mentioned m items 2,3, and 4, but wil 
mquimXque sensor on the second joint. The following analysis demonstrates 

the latter observation. 
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Figure 7.1 Proposed MPMS. 


Evaluate (5-2) and (5-3) for 04 - 90°, 

T keaI “ V c 2(01 C 2 - 2 0, e 2 s 2 ) + i„ 0 s 2 (2 0! 0 2 C 2 + 01 s 2 ) 

+ Ivzu (-02 c 2 + 02 2 S 2> - hm (®2 2 c 2 + 02 s 2> - hyu < 2 »1 C 202 + ®1 ^X 7 ' 1 ) 


\e2 = h,n 02 + ^yu Q l ' £ c 20, “ *zzu ®1 S 2 ^zu 6 1 C 2 


#i 2 s 20 2 + 5 V ®i 2 s a>2 (7 ' 2) 

Observe that I zzu cannot be determined from (7-1). However, if one measures x ke2 
and fixes 0^ =0, 0j = 0, then (7-2) yields 
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I«z. = ^ when 0,= 0,0, =0. 
0 2 


(7-3) 


To measure the other inertia parameters set 0 2 = 0, ®2 _ 0 in ( 7 "0 and ( 7 " 2 ^ to 
obtain 


U keul 


• ( lyyu C 2 


+ 1 


xxu 


S 2 2 " lx 


yu 


S20, ) 0 1 


(7-4) 


It I T 

(1-xyu C 29-> — 2 lixu S 20 7 ' 2 


yyu 


s 20 2 ) e i 


2 _ 


T ke2 = ( I xzu s 2 + *yzu c 2> ® 1 (7 ' 5 ^ 


The procedure would be to use (7-4) to obtain I yyu , Ixxu Jxyu an< ^ ^ en substitute 
these values in (7-5) to obtain the remaining inertia terms >Iy ZU • 


Experiments 7.1 and 7.2 entail measuring both torques, 0 X and 0 X with the second 
joint axis fixed, respectively, at 0 2 = 0, ^/ 2 and For ttie last angle ’ 
measurement of the torque on the second joint is not required. 


Experiment 7.1 

1. 0 2 = O , 


02 = 0,02 = ° 


Iyyu 


kel 


0 


2. 0 2 — Jt/2 , ^xxu 


3. 0 2 = jc/ 4 , I xyu 


'-kel 


©1 


^yyu C 2 *** ^XXU S 2 


- X 


kel 


9, 


(7-6a) 


(7-6b) 


(7-6c) 
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Experiment 7.2 



1. 0 2 - 0 , I yzu 




ke2 



(7-6d) 


2. 0 2 — tt/2 , I*^ 




ke2 



(7-6e) 


Equations (7-3) and (7-6) determines the inertia tensor of the unknown mass . 
Further, the same measurements performed on the MPMS when it is unloaded 
will allow one to compute the inertia tensor for the turntable and link. 


8. Conclusion 

Theoretical analysis proves the inertia parameters of an unknown mass can be 
determined from the joint torques, positions, angular velocities and angular 
accelerations of the Mass Properties Measurement System (MPMS). In particular, 
the existing system requires measurement of the torque and angular acceleration 
on joint one and the angular position, velocity and acceleration on joint two. A 
proposed system, where the twist angle between the MPMS turning axes equals 
90°, permits computing the inertia terms in a simple manner with the additional 
requirement of measuring the second joint torque for one set of measurements. 
The new design offers some significant advantages. It packages more compactly, 
allows simple mechanical scaling, and is easier and less costly to construct. 

The next report will discuss actual physical measurements with the existing 
system to determine the precision with which the MPMS can measure the mass 
properties of an unknown mass. 
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